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Abstract. — We revisit a long standing issue in the theory of disordered electron systems 
and their effective description by a non-linear sigma model: the hyperbolic Hubbard- 
Stratonovich (HS) transformation in the bosonic sector. For time-reversal invariant sys- 
tems without spin this sector is known to have a non-compact orthogonal symmetry Op^q. 
There exists an old proposal by Pruisken and Schafer how to perform the HS transforma- 
tion in an Op_ ^-invariant way. Giving a precise formulation of this proposal we show that 
the HS integral is a sign-alternating sum of integrals over disjoint domains. 



1. Introduction 

Initiated by work of Wegner [T|, Schafer & Wegner and Pruisken & Schafer ["35, 
non-compact non-linear sigma models as well as their supersymmetric generalizations 
due to Efetov [4J have long been a standard tool in the field of disordered electron 
physics. Areas of application include single electron motion in disordered and chaotic 
mesoscopic systems (jSJ, chaotic scattering |16li7l, localization and delocalization in 
systems exhibiting the Integer Quantum Hall Effect [|8l, statistical properties of the 
Dirac spectrum in non-abelian gauge field backgrounds [l9l, to mention only a few. 

The utility of the non-linear sigma model derives from the fact that it exposes the 
long-wavelength degrees of freedom of the disordered system, which are hidden in the 
original microscopic formulation by a Hamiltonian with quenched, i.e., time-indepen- 
dent, random parameters. These degrees of freedom correspond to interacting diffusion 
modes, known to be responsible for the universal behavior of spectral and eigenfunc- 
tion statistics in a broad class of disordered systems [HllSl. 

While the final results are expected to exhibit a high degree of universality, the math- 
ematical tools employed to derive the non-linear sigma model may vary depending on 
the type of microscopic model under consideration. A concise introduction to this issue 
can be found in [|TOl ; see also the recent papers lfIT| [T2l. For some types of microscopic 
model the following mathematical identity plays a key role in the derivation: 
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where Cq is a constant and \dR\ denotes the Lebesgue measure of a normed vector 
space of matrices i? . In a situation with compact orthogonal symmetry, which arises 
when the fermionic sector of the theory is considered, one has license to take the in- 
tegration domain D to be the real symmetric matrices R, and it is then trivial to do 
the Gaussian integral by completing the square and shifting variables. However, for 
reasons reviewed in Sect. [2] the bosonic sector (say, of time-reversal invariant systems 
without spin) calls for A to be a non-symmetric matrix composed of elements 

^U=l!l=ina(Pj,aSj, (1.2) 

with (p + q) xN real numbers (pi, a and si = . . . = Sp = — 5p+i = . . . = —Sp+q = 1 . 

The diagonal matrix s with entries 5// = Si 5ij determines a non-compact variant ^ 
of the real orthogonal group by the condition g^sg = s , where g^ means the transpose 
of the matrix g. First discovered in the present context by Wegner the group Op^q 
is sometimes referred to as the 'hyperbolic symmetry' of the problem at hand. 

In the situation with hyperbolic symmetry, choosing a good domain of integration 
fori? and making rigorous sense of the integral (11.11 ) are non-trivial tasks. The difficulty 
arises from the general context of the so-called Hubbard-Stratonovich method: the 
exponential e^'^'"'^^ must be kept bounded - or else the next step of the method, 
which is to integrate over the microscopic fields , would be invalid. Since the real 
matrix A satisfies the symmetry relation A = sA^s one might think that one should try 
to integrate over the domain of all real matrices R subject to the symmetry R = sR^s. 
Unfortunately, the resulting quadratic form TrR^ = Tr RsR^s is of indefinite sign and 
therefore such a choice of integration domain for R makes the i?-integral divergent. 

Nonetheless, a valid solution of the problem posed by (11.11) . i.e., a choice of integra- 
tion domain for R making the 7?-integral converge while keeping the integrand bounded 
as a function of A, was offered in the paper by Schafer and Wegner [|2||; by deforming 
some of the real freedoms in R into the complex numbers they parameterize R as 

R = P-i?iTsT^ , 

where P runs through the real symmetric matrices that commute with s , the constant A 
is any positive real number, and T G SOp q. The reader is referred to the review [HOl 
for a detailed discussion. Now the Schafer- Wegner choice of domain has, besides its 
many merits, a certain drawback: it lacks invariance under the action of the symmetry 
group Op,q. This may be the reason why the Schafer- Wegner choice was never much 
in use by the disordered electron physics community. 

1.1. Statement of result. — Another choice of integration domain for (11.11) . which 
meets the afore-mentioned requirements and also has the desirable property of Op^q- 
invariance, first appeared in a paper by Pruisken and Schafer [3J. Beginning with ^ 
this choice came to be used in numerous later works. In the present paper we are going 
to establish this so-called Pruisken-Schafer domain as a rigorous alternative to that of 
Schafer- Wegner. There exists, however, a mathematical subtlety: although our choice 
of domain agrees with Pruisken-Schafer as a union D = Ucr^cr of sets Da , it differs as 
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an element = Ecr sgn(a)Dcj of the vector space of integration chains (for the details, 
see below). Some low-dimensional special cases of it were worked out in [[I3t [l4 l| . 

The Pruisken-Schafer domain, D, is an open subspace of the j{p + q){p + q + l)- 
dimensional space of real matrices R obeying the mixed symmetry relation R = sR^s. 
Because such matrices may in general have complex eigenvalues and eigenvectors, one 
defines D as the subspace of matrices R = sR^s that can be diagonalized by conjugation 
R g^Rg with an element g E Op^q of the real (non-compact) orthogonal group. 
This set D turns out to be the union D = \J^Da of {^p) = (^+^) sub-domains Da . 

The domains Da are enumerated as follows. With the exception of cases forming a 
set of measure zero, an Op^^-diagonalizable matrix R = sR^s has p 'space-like' and q 
'time-like' eigenvalues A G M - here we use the language of relativity theory to com- 
municate that the corresponding eigenvector v has the property v'^'v > or v'^'v < 0, 
respectively. In the generic situation without degeneracies, the eigenvalues of R can 
be arranged in decreasing order and this ordered sequence translates into a binary se- 
quence o by writing, say, the symbols • for space-like and o for time-like eigenvalues. 

The binary sequence o encoding the relative order of space-like and time-like eigen- 
values is an Op, ^-invariant. Moreover, collisions between eigenvalues of space-like and 
time-like type generically lead to the birth of complex eigenvalues, thereby taking us 
outside the integration domain for R . This means that the Op.^-diagonalizable matrices 
R of binary sequence o form an Op, ^-invariant domain, which we denote by Da , and 
any two such domains Da and Da' with o ^ o' can only touch each other in subspaces 
of lower dimension (more precisely, of codimension two). 

The integration measure \dR\ is defined to be the fiat one for all domains Da '■ 

\dR\=ll^^jdR,j. (1.3) 

Now the startling feature of the following statement is that the integral on the right- 
hand side of (11.11 ) is proposed to be a sum of integrals /^^ with alternating sign! 

Theorem 1. — There exists some choice of cutojf function R ^— Xe{R) (converging 
pointwise to unity as £ ^ 0), and a unique choice of sign function o sgn(c7) G {±1 } 
and a constant Cp_q such that 

Cp,, lim£sgn(a) / ^-^'^"-^^'^^Xe{R)\dR\= (1.4) 

holds true for all matrices A = sA^s with the positivity property As > 0. 

Remark. — It will be shown that sgn(a) is the parity of the number of transpositions 
• ^ o needed to reduce the binary sequence a to the extremal form Oq = • ■ ■ ■ ■ ■ ■ o. 
It should be emphasized that (11.41) is not a standard Gaussian integral; in fact, we have 
not succeeded in finding a proof of this formula by completing the square and shifting. 

It is informative to give an alternative expression for the integral (|1.4I) in terms of the 
eigenvalues of R . Let R = g?ig^^ with A = diag(Ai , . . . , ^p+q) G RP+^ and g G SOp,^ . 
The volume element transforms as \dR\ = J{X)\dX. \ dg where \dX \ = Y[i=\ dXi and dg 
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is a positive Haar measure for the connected group SO^^^ . The Jacobian is 

Corollary 1. — Let a function /'(A) with alternating sign on W^^^ be defined by 

/(A)=7(A)n n sign(A,- - A,) . 

Then, assuming the positivity As > of the matrix A = sA^s, we have 

Cp,,lim / ([ c-^-^^^^^'^-'xeis^g-')dg)c-^^''j'iX)\dX\=c-^^^\ 

Remark. — The first p eigenvalues Xj are space-like, the last q time-like. To deduce 
Cor. [Hfrom Thm. [H we need only observe that the function J'{X) differs from the 
Jacobian 7(A) precisely by our sign function A t-^ sgn(c7(A)): 

/(A)=7(A)sgn(a(A)). 

The presence of this alternating sign factor was conjectured in [|l4ll : see also lll3ll . 

In previous work, beginning with [6J, the identity of Cor. [T] was assumed to hold 
without the alternating sign. It became apparent in 1995 that this assumption was un- 
founded [15] though it remained unclear at that time how to correct the mistake. Prob- 
lems with the Pruisken-Schafer domain were mentioned in lfl6l . but for various reasons 
the issue was never much emphasized in the published literature. (For one thing, all 
results derived using y(A) instead of /(A) stand correct after saddle-point approxima- 
tion for large A^, as the integral in the large-A^ limit is dominated by contributions from 
a single domain Da ■ For another, the Schafer-Wegner domain was available [16] as a 
rigorous alternative in order to do all intermediate steps correctly.) The present paper 
solves this long standing problem in a satisfactory if perhaps surprising manner. 

Let it be pointed out that the matrix As appearing in applications is not strictly posi- 
tive but positive 5em/-definite. This is, however, a minor issue as an easy variant of the 
formula (11.41) takes care of the semi-definite case; see [[l3ll for the details. 

For the related case of non-compact unitary symmetry Up,g a formula analogous to 
that of Cor. [H had been established in [13], using the special feature of semiclassical 
exactness by the Duistermaat-Heckman theorem [17]. Although it will not be shown 
here, the methods of the present paper are robust and can be adapted to handle the case 
of Up^^ as well, without taking recourse to semiclassical exactness. 

The paper is organized as follows. In Sect. [2] we present some background mate- 
rial concerning the Hubbard-Stratonovich method. Basic results needed from integral 
calculus are collected in Sect. [3l Then, in Sect. [4l we work out the simple but instruc- 
tive case of p = q = I in detail. The case of general p and q is handled in Sect. [5] by 
reduction to p = ^ = 1 . 
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2. Background 

To make the present paper self-contained, we now describe the steps prior to (11.41 ). 
called the Hubbard-Stratonovich transformation in the present context. Although the 
true power of this transformation is to allow the treatment of non-trivial models with 
local Om gauge symmetry [12|| or even without such symmetry [3], we will restrict 
ourselves here to reviewing the basic steps at the example of the simplest random 
matrix model of this universality class: the Gaussian Orthogonal Ensemble (GOE). 

In the present work we are concerned with time-reversal invariant systems without 
spin. For such systems, the quantum mechanics of stationary states can be done over 
the field M of real numbers. Assuming the real Hilbert space to be of finite dimension, 
A^, we may express the Hamiltonian operators H by real symmetric N x N matrices. 
The Gaussian Orthogonal Ensemble by definition is an OA^-invariant probability mea- 
sure dil{H) on the space of such matrices. For our purposes, the GOE measure dil{H) 
is best characterized by its Fourier transform: 

(^'^''''')goe ^ = e"^^'"''' ' (2-1) 

where the Fourier variable K is any symmetric matrix, and G M is a parameter. 

An important object of the theory (see, e.g., [lOJ) is the expectation value of the 
reciprocal of a product of (square roots of) characteristic polynomials of H : 

F(zi,...,z,+,) := (n;;f Det-^/^Zj-H))^^^ (z, G C\M) . (2.2) 

We assume that 3mzj > for j = 1 , . . . , p and '^mzj < for j = p + I, . . . ,p + q . 
To compute such an expectation value, we use the trick of representing the reciprocal 
square root of each determinant as a Gaussian integral over vector variables (p : 

where \d(p\ = Y[a=i^9a and the factor sj := sign(Jmzj) ensures the convergence of 
the integral. '■= La=i V^a means the Euclidean scalar product of M^. 

Substituting such integrals into (12.21) . one immediately performs the GOE average by 
employing the characteristic property (12.11) . The result can then be expressed in terms 
of a matrix A{(p) of size {p + q) x {p + q) with entries A{^)ij = Y.a=i ^',a9j,aSj 

F{zu...,z,,,)= /e'^P?'"-^<^'-''''e-|J™(»)=n7rJ^. (23) 

J j [iTtSjj 

Here is where the identity (11.11) comes in: it enables us to linearize the quadratic term 
Tr A{(p)^ in the exponent, thereby reducing the integral over the A'^-component vectors 
9i , . . . , cpp+q to decoupled Gaussian integrals, one for each component a = I,. ..,N. 
When the latter integrals are carried out, (12.31) becomes an integral over the collective 
variables R only. Using the rigorous form (11.41 ) of the relation (11.11) we obtain 

F{z)=C'p^,Y.^gn{a) f Q-^^''''^a-'^'\z-R)\dR\, (2.4) 
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where z '■= diag(zi, . . . and, with the assumption that A^^ is large enough, the 

cutoff function Xe has been removed by sending the regularization parameter e to zero. 

The integral representation (12.41) for F{z) is well suited for saddle-point analysis in 
the limit of large A^. We can now appreciate why the identity (11.11) is useful here: it 
serves the purpose of exposing the good variables in which to perform the large-A'^ 
limit. Let this fact suffice as a motivation for our labors in the body of this paper. 



= . (3.2) 



3. Basics from calculus 

Recalling from Sect. II. li the definition of the connected domains Da (cf. also below) 
we consider the following alternating sum of integrals: 

/(A)= lim ysgn(cT) / t-^'^''+''^^\e{R + -iA)\dR\ , (3.1) 

where is some smooth cutoff function which regularizes the integral and converges 
pointwise to unity in the limit e — > . The real matrix A is subject to the conditions 
A = sA^s and As>Q. Choosing the Lebesgue measure \dR\ as in (|1.3I) we shall prove 
that the integral /(A) does not depend on the matrix A, or equivalently, that all of its 
directional derivatives vanish: 

Once (|3.2I) has been established, the Hubbard-Stratonovich transformation (|1.4I) fol- 
lows by multiplying both sides of (13.11) with a constant e^^'^'^ . 

To prepare our treatment, we recall a few basic facts from calculus. Given a vector 
field V on an n-dimensional differentiable manifold M, let [—5,5] 3 t ^ be the flow 
of V, i.e., the one-parameter family of mappings (pt'.M^M determined by ^f=o(i?) = P 
and -^^tip) = v{^tip)) for all p e M.In local coordinates x^,...,x",we have 

^x'o^t=v'o^t {v = v'd/dx') . 

We are using the summation convention: an index appearing twice (once as a covariant 
and once as a contravariant index) is understood to be summed over. 
For a differential form a = oc,,.,.,^^ dx'^ A ■ ■ ■ A Jx'* on M let 

= (a,i.../i o ^f) d{x'^ o^f) ■ A o ^f) 

denote the puUback of oc by ^/ . The Lie derivative of a w.r.t. the vector field v is then 
defined by differentiation at f = : 

dt /=o 

By the relation A /3) = a) A {^* {5) , the Lie derivative satisfies the Leibniz 
rule ^v(a A /3) = {^vU) A/3 + a A (=5fv/3). Cartan's formula expresses it as 
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where d is the exterior derivative, 

d (a,j.../^<iyi A- - ■ A jy*) = dai^,„i^Adx'^ A ■• • A Jx'* 
and ly denotes the operation of contraction with the vector field v : 



(iva)/,.../,_, = V («//,.../,_, - a,i/...,,_; + . . . + (-1) 



After this summary of basic facts, let to be an n-form on M and consider the integral 
jf^{^vf)o} where = i^df is the derivative of some function / : M ^ C in the 
direction of the vector field v. Assuming that a is invariant by the flow of v, the 
Leibniz rule gives 

JM Jm Jm 

By Stokes' theorem for the integral of the exact n-form d{flvCo) it follows that 



(^v/)«= / /iv«, (3.3) 

M JdM 

where dM is the {n — 1) -dimensional boundary of M. 

In what follows, we will employ (|3.3I) in the closely related case where one is inte- 
grating against a density \o)\ instead of an n-form CO. Formula (13.31) still holds true in 
that case - assuming of course that Jfy\0}\ = 0. Note also that in order to integrate the 
twisted {n — l)-form fiy\(o\, the boundary dM must be equipped with an outer orien- 
tation, i.e., a transverse vector field pointing, say, from the inside to the outside of dM. 
(Technically speaking, a density |ft)| on M is a top-degree differential form co tensored 
by a section s of the orientation line bundle of M, and a choice of splitting \ (o\ = CO0s 
determines an isomorphism between outer and inner orientations for dM.) 

To apply the formula (13.31) to the situation at hand, let t{R) = Rjj d / dRij denote the 
vector field generating translations 0f (i?) =R + tR. Then 

^HA + tA) = lim i£sgn(c7) / ^^^Jc-^'^''^^^' MR + iA)) \dR\ , 



where t(A) = Ajj d / dRij . Here, recognizing the fact that our integrand depends on R 
and A only through the combination 7? + iA , we have used the identity 

d ... . d 



/(i? + iA)=i3— /(i? + iA). 



<9A,/ ' ' dRij 



Now the Lebesgue measure (or positive density) \dR\ is invariant under translations 
R^ R + tA. Therefore |^^^| =0, and we may apply formula (13.31) to obtain 



d 

dt' 



I{A + tA) = lim iV sgn(a) / Q-^'^''+''^'^\e{RViA)i^a)\dR\ . (3.4) 



The main achievement of this paper is a proof that the alternating sum of integrals 
on the right-hand side vanishes in the limit e — > . For that purpose we need to develop 
a good understanding of the boundary components dDa and the twisted form i^f^j^^-^ \dR\ 
restricted to them. As a first step, we illustrate the essential features of the argument at 
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the example of the symmetry group being Oi i . This example plays a very important 
role as the general case will be handled by reduction to it. 



4. The case of O i i -symmetry 

Here we are going to deal with the special case of real 2x2 matrices R subject to 
the linear symmetry relation 

R = sR^s, 5 = diag(l,-l) , 

where R\-^ R^ means the operation of taking the matrix transpose. Such matrices can 
be parameterized by three real variables rn , ri2 and as 



R 



m ri2 
-rn r22 



Thus our space of matrices R is isomorphic as a vector space to ^ . The measure of 
integration is the positive density \dR\ = dr\\dr\2dr22 ■ The matrix A is of the same 
form as R but its matrix elements are constrained by 



an > > a22 , \an\ < V-ana22 



ail <^12 
-ai2 a22 



Another way of describing the setting is to equip the real vector space with the 
sign-indefinite bilinear form B determined by s : 

B{u, v) = u^sv := uivi — U2V2 ■ 

The matrices R are then symmetric with respect to B in the sense that 

B{u,Rv) =B{Ru,v) 

for all a,v G M^. We denote the linear space of real matrices R with this symmetry 
property by Symg(]R^). Since B has signature (1,1), the symmetry group of B is the 
non-compact real orthogonal group Oi 1 . Elements ^ G Oi 1 satisfy the equation 

s = g^sg , 

which is equivalent to saying that v) = B{gu,gv) for all w, v G M^. 

For reasons indicated in the Introduction (see Sect. [T]) - let us recall that in order for 
the integral (|l.ll) to exist one needs the inequality TrR^ > , at least asymptotically - 
we want our matrices R G Symg(M^) to be diagonalizable by the action R t-^ gRg^^ 
of the real group Oi 1 by conjugation. The condition of diagonalizability is formulated 
most clearly by expressing R 'ma basis of light- like vectors of M?: 



e+ :-- 



0' ^-^=(-1 



which obey 5(e+ , 6+) = B{e , e-) =0 and B{e-^ , e-) = 2. Applying R to this basis 
one has 



RIGOROUS HUBBARD-STRATONOVICH 



9 



where A = 5(^11 +r22), r] = 5(^11 -r22) +^2, and ^ = \{rn-r22) -rn - Since the 
diagonal piece of R in this basis is a multiple of unity, the condition for diagonalizabil- 
ity is that the product of off-diagonal matrix elements be positive: 77 > . 

The domain ^r] > of Oi i-diagonalizability of our matrices R E Symg(M?) consists 
of two connected components: 

D.O : <^,T] > 0, Do. : <^,T7 < 0. 

The motivation for the notation (•, o) will become clear shortly. Here we simply remark 
that if V,, Vo denotes a pair of ^^-eigenvectors with 5 (v., v.) > > 5(vo, Vo) and A,, Ao 
are the associated eigenvalues, then A, > Ao for R e D,q and Ao > A, for R e Do. . 
Next, let Ga^e denote the regularized Gaussian integrand 

G^,,(i?):=e-T'-(^+^)^;^,(i? + iA), 
where we choose the cutoff function to be 

;t,(7?):=e-iTr(^-^-^.0^=e-4^^'->2'. 
Using the coordinates A , , 77 we have the expression Ga.e = e^-^2-/i -/o ^jt^ 

/o = +ay +2^22(1 -2e) , /2 = 2{X^ + ^ri) +e{ri -^f , 
/i =2i(aii-fa22)A-fi(aii-a22)(^ + ^7) + 2iai2(l-2e)(^-7]) . 
Let us observe that the integral /(A) defined in (13.11) can now be written as 

4.1. Description of boundaries. — By their definition via inequalities, the domains 
Da (for o = •o,o«) are open in the three-dimensional space Symg(R2) and have 
boundaries dDa which are two-dimensional. Each of the boundaries dDa is a union 
of 2 half-planes. To describe this in detail, let 4 half -planes Ea and Fa be defined by 

£.0 : <^ = , T] > , F.o : T] = , ^ > , 
£0. : <^ = , 77 < , Fo. : 7] = , <^ < . 

Ea and Fa make up the boundary of Da (for both o = •o and o = o»). Being part 
of dDo» , the half -plane Eo» inherits from Do, an outer orientation by the transverse 
vector field = d /d^ . (Indeed, starting from any point R E Do* very close to the 
boundary Eo, and making a small step in the direction of , one crosses Eo» .) In the 
same sense, Fo, C dDo, is oriented by drj while Ea C dDa and Fa C dDa for C7 = #0 
are oriented by the opposite vector fields —d^ and —dri respectively; see Fig.[Tl 
With these orientations being understood, we have the relations 

dDa=Ea+Fa (a = «o,o«), 

which are equalities in the sense of chains, or distributions. (To achieve equality in the 
sense of sets, we would have to include the real line = rj =0, which is part of both 
(9D,o and ^Do, . However, our aim is to integrate bounded differential forms, and for 
that purpose the lower-dimensional parts of the boundary are of no concern.) 
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Figure 1 . Each sub-domain Da has two boundary components, E(j and . 
The outer orientation of these boundaries is determined by the rule of crossing 
the boundary "from the inside to the outside". 



With dD(y as described above, formula (13.41) now reads 

^liA + tA) = \im ( f -f ]iGA,ei,(A)\dR\- (4.1) 

4.2. Reorganization of boundary pieces. — The key idea is to reorganize the bound- 
ary half-planes into consistently oriented closed manifolds. Define the full planes 

E : ^ = (oriented by —d^), 
F : rj =0 (oriented by — drj) , 

so that we have 

E = E^o Eo» , F = Ffo 7 

in the sense of chains. Since 

dD,o - dDo, = E,o + F.o -Eo,-Fo,=E + F 

(still as chains), we obtain the following relation between integrals: 



CO- C0= C0+ CO 

dD.o JdDo. Je Jf 

for any twisted 2-form CO which is bounded (if not continuous or smooth). In particular, 
this relation holds for co = iGa^e It(a) I^^I- Thus (14.11) turns into 
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F 



Figure 2. Reorganization of the boundary pieces: defining E = E,o —Eq 
and F = F,o — Fo, one obtains two consistently oriented 2-planes E and F. 



A graphical sketch of the situation is shown in Fig. [21 

4.3. Computation of boundary integral. — We are now going to show that the 
integral co vanishes for O) = iGa^e ^^(a) I^^I '^he limit e ^ . Essentially the same 
calculation shows that jp (o vanishes in this limit. 

To handle the case of let i// : £ ^ Sym5(M^) be the identical embedding and 
introduce coordinates A , T] on £ so that 

Recalling Ga.e = e^^2-/i-/o back the functions /2 and fi to E: 

\j/*f2 = 2X^ + £ri^ , v/7i =2i(aii+a22)A+i(flii-a22-2ai2(l-2e))77 . 

Next we compute the puUback of the twisted 2-form i^^^^ \dR\. For this, we write 

\dR\ = dr\idr\2dr22 = dXdt,dr\ := {dX Ad£, A Jt]) ® Or , 

where Or = +1 when Symg(M^) is oriented by the ordered set of linear coordinates A, 
r] (or any even permutation thereof), and Or = — 1 when the ordering X, rj, ^ (or 
any other odd permutation of A, i^, r]) is chosen. Using the coordinate expression of 
the vector field t(A), 

T[A) = an^ hai2-5 l-a22 



drii dri2 dr22 ' 
\{dii+d22)dx + \{dn - d22 -2di2)d^ + j{dn -^^22 + 2^12)^7] 
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we then obtain 

V/*i^(^)|Ji?| = j{du -^22 -2ai2) (Jt] AJA) (g)Or . 

Now the presence of the orientation factor Or tells us that an outer orientation by ±d^ 
translates into an inner orientation by ±Jt] AdX. Since our At] -plane E is oriented by 
—d^ , we must assign to it the orientation given by —drj AdX . In this way, assembling 
everything we arrive at the expression 

/ rGA,eh(A)\dR\ = -lian -d22 -2di2)e«n+"22-2«?,(l-2£) ^ 
•J E 

Jr Jm 
where the right-hand side has been reduced to a product of Riemann integrals by Fu- 
bini's theorem and the definition of what it means to integrate a differential form. 

The crucial point now is that the integral over rj vanishes in the limit e — > . Indeed, 
writing be := ^(aii —6(22) — 6(12(1 — 2e) for short, the T]-integral gives 



f e-^'j'-2i*^Vr]=e-*'/^v^ 
Jr 



and this does go to zero for £ ^ as long as the real number bo is non-zero. The latter 
condition is always satisfied, as the constraints on the matrix elements of A ensure that 

bo = j{ai\ -aji) -an > V-anoxz- {anl > . 

Thus we have shown that lime^o+ Ie ¥*^Ga,£ It:{a) l^^l = ^ • 

The situation is no different for the integral over the 2-plane F. We therefore con- 
clude that ^I{A + tA) |/=o = . Hence /(A) =: Cj^J is a constant independent of A for 
As > 0. This proves the validity of the Hubbard-Stratonovich transformation in the 
following precise form with cutoff function Xei^) = e^i'^'^i^^^^^') : 

e-Ti-^'=Ci,i lim y sgn((7)/ e-^''''-^'^'^^XemdR\ , 

where sgn(«o) = 1 and sgn(o«) = — 1, and we have replaced Xe{R + iA) by Xe{R). 

Let us emphasize once again that the right-hand side is not a standard Gaussian 
integral, and that our proof of this formula was not by completing the square and 
shifting. We will show in the next subsection that C14 = 122 ;r 2 . 

4.4. Short proof. — We now spell out another line of approach for the case of ;? = 
q= I, which is a variant of that given in llT3l . This calculation will be much quicker, 
but does not extend (at least not in any way known to us) to higher values of p, q. 
Recall that in light-cone coordinates ^, rj, and A, we have the expression 



where bo = j{aii —aji) — ai2 > and bi = j{a\i — a22) + «12 > 0. It is clear that this 
can be integrated against dX resulting in a factor a/tt e~ ^ ("1 1 +"22)^ \[ tj^g^ remains to 
integrate the other factor Q^^''^^^^ori--^b\ t, against dt, dr\ over the two quadrants , r] > 
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and (^,7] < 0. For this we make the substitution of variables r/ = pe^ and ^ = pe^ 
which transforms the measure dr\ to 2|p | dpdx . Noticing that p > on D,o , p < 
on Do, , and our chain of integration is D = Z),o — Do* , we drop the absolute value on 
IpI and integrate against the sign- alternating distribution /]j2 pdpdx . More precisely, 
we take the p integration to be the inner one (this obviates the need for regularization 
albeit at the expense of Fubini's theorem becoming inapplicable, so that the order of 
integrations is now fixed and can no longer be changed). The resulting p integral is 

Jr 

We now observe the relation 

/3e-^'jT = e-*"^i e-5(^oe^-*ie-^)2 1 J (^Qgt ^-t^ ^ 

which shows that the remaining T integral is another standard Gaussian integral in the 
variable j{boe^ — hi e^^), yielding the value v^e^*"^' . Altogether we then obtain 

' Jd 

for D = D,Q —Do,, which is the desired result but for scaling A —>■ \f7.A and R \flR . 
(Of course it must be understood here that the p integration has to be done first.) 

This computation, while direct and short, does not seem to carry over to the case of 
higher values of p and q, whereas our conceptual proof of Sects. l4~Tll4.3l does. 

5. The general case of 0^, ^-symmetry 

We now turn to the general case of real matrices R of size (p + q) x (p + q), 

R = ( ^P^P ^P'l\ 
\-{Rp,qy Rq,q) 

where the blocks Rp^p and Rq^q are symmetric. Such matrices i? as a whole satisfy the 
linear symmetry relation 

R = sR^s , s = diag(Idp , — Id^) . 

Equivalently, if we equip the real vector space W^'^ with the indefinite bilinear form 

p p+1 

1=1 j=p+i 

then our matrices R form the linear space Symg(M^+'5') of 5-symmetric matrices, i.e., 

Vm,vgM^'+^: B{u,Rv) = B{v,Ru) . 

The symmetry group of the bilinear form B with signature {p,q) is the non-compact 
real orthogonal group O^^^ of matrices g satisfying the equation s = g^sg or equiva- 
lently, 5(m,v) = 5(^i<,^v)'for all M, V G M^+^. 

We still want our matrices R G Symg(R^'+^) to be diagonalizable by the O^.^-action 
R i-^ gRg^ . In the simple case p = ^ = 1 we saw that diagonalizability fails when 



14 



Y.V. FYODOROV, Y. WEI & M.R. ZIRNBAUER 



the off-diagonal block Rp q is larger than the difference of the diagonal blocks Rp^p 
and Rq.q- A similar phenomenon is expected to occur for general p,q. Roughly 
speaking, R will fail to have real eigenvectors and eigenvalues when the degree of 
non-Hermiticity coming from the off-diagonal blocks becomes too large. 

5.1. Description of components D(, . — In the Hubbard-Stratonovich integral (13.11) 
we integrate over the set, D,of R G Symg(]R^+^) that can be brought to diagonal form 
by the Op,^-action, i.e., matrices of the form R = gX with A = diag(Ai, A2, . . . , ^p+q) 
and g E Op^q. We know from our considerations for p = ^ = 1 that D decomposes into 
connected sets Da whose boundaries intersect in lower-dimensional sets (of codimen- 
sion 2). Here we describe these connected components Da ■ 

Let R — gXg^^ be Op ^-diagonalizable in the sense just described, so that the p + q 
column vectors which make up g form a basis of ]RP+^ consisting of eigenvectors of 
the matrix R . Let v*^^) denote the column of g. Then by the equation g^sg = s the 
vector v*^^) has positive norm square B{v^^\v^'^^) = 1 >Oforfc = l,...,p and negative 
norm square B{v^^\v^^^) = — 1 < for k = p + I, . . . ,p + q . To distinguish between 
these two cases, we use the language of relativity theory and call the corresponding 
eigenvalue Xk of 'space-like' type in the former case and 'time-like' in the latter. 

Assuming that we are in the generic situation, where all eigenvalues Xk of R differ 
from each other, we can arrange them as a descending sequence. Given any such 
sequence, we assign to it a motif, a: if an eigenvalue of the ordered sequence is space- 
like, then we represent it as a full dot, otherwise as an empty dot. Thus our motifs o 
are sequences of p full and q empty dots; e.g., c7 = «oo«o«isa motif for p = q = 3. 

With each motif o we now associate a domain Da as follows. While insisting that 
eigenvalues of opposite type (space-like vs. time-like) must never be equal, we lift the 
condition of non-degeneracy of eigenvalues partially: eigenvalues of the same type are 
now allowed to collide and exchange positions, but eigenvalues of opposite type still 
are not. This rule of eigenvalue motion defines a connected domain Da for each a . 

We define a sign function a t— > sgn(a) G {±1} by setting sgn(c7) := (—1)" where n 
is the number of transpositions •© ^ o» which are needed to reduce o to the reference 
motif C7o = ««---«oo---o made of q empty dots following p full dots. The parity of 
this number n is well-defined although n itself is not. Indeed, another expression for 
sgn(a) = (-1)" is 

p p+i 

sgn(c7) = n n sign(A/ - A/) , 
where the A 's are the eigenvalues of any R G Da . 

5.2. Description of boundary. — As was explained above, a point in the interior 
of Da is a matrix R E Symg(R^+^) which is diagonalizable by the Op^ ^-action and 
has eigenvalues that order according to the motif C7 - in particular, no two eigenvalues 
of opposite type are equal. This means that if t h-^ R{t) is a continuous curve which 
starts in Da but leads to a collision of two or more eigenvalues of opposite type at 
tc, then R{tc) is a point in the boundary of Da ■ The generic boundary situation is 
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that (i) exactly two eigenvalues of opposite type collide, (ii) the norms of the two 
corresponding eigenvectors go to zero, and (iii) the boundary point R{tc) is a matrix 
which fails to be diagonalizable by any real transformation g E Op q. 

To carry out the computations leading to the desired result (13.21) . we need a good 
description of the interior of D = IJa ^cr as well as the boundary dD. This is achieved 
as follows. Let us first make a slight change of perspective: abandoning the viewpoint 
of i? as a matrix, we switch to regarding R as an invariantly defined linear operator 
(giving rise to a matrix when expressed w.r.t. some fixed basis); i.e., taking the matrix 
elements R'j we let 7? := R'-^Ejj (summation convention!) with Ejj the standard basis 
of End(R^+^) defined in terms of the standard basis e/ of RP+^ by Eij = e,- dji^ . 

Consider now the Grassmannian Gri i(M^+^) of 5-orthogonal decompositions 

where V C R^+^ is any Lorentzian plane, i.e., a 2-dimensional subspace V = L+ ©L^ 
spanned by a pair of null lines or light-like directions such that s{L^) = . The 
restriction of 5 to V is then non-degenerate, and the 5-orthogonal complement V"^ of 

V is spanned by p — 1 space-like and q—\ time-like vectors. Note that Gri i(R^+^) is 
acted upon transitively by the group SOp.g , which is to say that every Lorentzian plane 

V C R^'+^ can be regarded as the image ^(R^+^) of the standard Lorentzian plane R^+^ 
by some transformation g E SOp^q. Since R^+^ has isotropy S(Oi^i x Op-i^q-i), our 
manifold Gri i (R^^+^) is the base of a principal fibre bundle 

S0p,^^S0p,g/S(0i,i xOp-i,^-i) ~Gri,i(R^+^) . (5.1) 

Moreover, since the Lorentzian planes V C RP+^ are in one-to-one correspondence 
with pairs of real lines, one in R^^ and R^ each, Gri j (R^+^) is a product 

Gri,i(R^^+^)~RPp_ixRP^_i 

of real projective spaces RP,i-i = R"/ (R\ {0}) for n = p and n = q. From this identi- 
fication it follows immediately that the manifold Gri i(R^+^) is compact, non-simply 
connected, and non-orientable unless both p and q are even. It is also seen that (|5.1I) 
decomposes as a direct product of two principal fibre bundles 

S0„ SO„/S(Oi X 0„_i) RP„_i 

for n = p, q. We note that these principal fibre bundles are non-trivial. 

Choosing a point x E Gri i (RP+^) is equivalent to choosing an orthogonal projector 

By orthogonality, B restricts to a bilinear form of signature (1,1) on and a bilinear 
form of signature {p — l,q — I) on the orthogonal complement V^^ = (Id — n,^)RP+^. 
We denote by 0{Vx) resp. 0{V^) the symmetry groups of these symmetric bilinear 
forms. Let Symg(y,;) denote the subspace of linear operators r : Vx ^ Vx which are 
symmetric with respect to B restricted to V^, and let D{V^) C Sym^{Vj^) denote the 
subspace of 0(y^^) -diagonalizable linear transformations t : Vj^ Vj- . 
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Now consider the set, of triples (jc;r,f) consisting of any point ;c G Gri i(]R^+^) 
and two linear operators r E Symg(V'x) and t E D{Vj-). This set ^ has the structure of 
a fibre bundle 

;r: ^Gri,i(M^+«) , {x;r,t)^x, 
where the fibre over jc is a direct product 

One may also view ^ as an associated bundle = Gx^W with G = SOp x SO^ , 
K = S(Oi X Op_i) X S(Oi X 0^_i) , W = SymB(Mi+i) x D{R^P-^'^+^'i-^'^) , 

where the action of ^ on is by conjugation. Since the principal bundle G ^ G/K = 
Gri^i(M''+^) is non-trivial, so is the associated bundle 3^ Grij(MP+^). 

can be viewed as a subset of our space Symg(M^+^) of operators R . Indeed, for 
R = {x;r,t) we may use the projector 11;^ to define R v for any vector v E by 

Rv:=rnrV + t{ld-nj,)v . 

In other words, R acts on Vx as r and on as t . Since both r and t are symmetric 
with respect to B, so is i? = rn^ + ? (Id — IT;,;). More formally, this correspondence 
{x;r,t) rUx + tild-Tlx) defines a mapping ^ : 3^ ^ SymgiW^'i). Note that the 
mapping is not surjective, as the operators R with more than two non-real eigenvalues 
cannot be presented in this form. However, the image ^ ( ^) is large enough to contain 
our integration domain D and also (in the measure-theoretic sense) the boundary dD. 

Now for any point x of the Grassmann manifold Gri i (R^+^) let denote the two 
null lines of the Lorentzian vector space ■ In other words B{u, u) =0 for u E or 
u eL^ . The three-dimensional space Symg(yx) then decomposes as 

SymB(y,) = Mldy, ©Hom(L+,L-) ©Hom(L-,L+) . (5.2) 

Let us choose a basis ef E Lf of null vectors of so that se^ = ef and e^) = 1. 
In such a basis IT;,: is expressed as 11^ = e^B{e~ , ■) + ^x^i^t^ ')■ Linear coordinates A , 
^, rj for Symg(Vx) are then introduced by decomposing r according to (|5.2I) : 

r = A n, + ^ 4^(4, ■) + r? •) . 

Extending the notation used before, let Ex E Symg(V';i) and Fx C Symg(yx) denote 
the linear subspaces which are defined by the equations ^ = and f] =0 respectively. 
These 2-planes Ex and Fx will assume the roles played by E and F in Sect. 14.21 It should 
be emphasized, however, that ^ and rj do not extend to globally defined coordinate 
functions for 3^ . The point here is that transporting the pair of null lines Lf around 
a closed path in Gri i(M'^+^) one may find that L+ and get interchanged. In other 
words, although Lf are invariantly defined as a pair of lines, there is no invariant way 
of telling which line is L+ and which is . Asa result, the linear coordinate functions 
^ and rj projecting from Symg{Vx) to Hom(L^,L+) and Hom(L+,L^) are only defined 
locally in the variables E Gri i(M^'''^). By the same token, the union of planes ExUFx 
is invariantly defined for all x, whereas the planes Ex and Fx individually are not. 
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The smooth assignment Symg{Vx) D E^UF^ x gives a bundle S' C ^ with fibre 
{E^UF^)xD(y^^) overxGGri,i(MP+^). By applying the map : ^ Symg(MP+^) 
we then obtain a submanifold (j){(S') C dD. As a matter of fact, since the generic 
situation at the boundary of D = D(M^+^) is that exactly two eigenvalues of opposite 
type collide, the set ^(co') agrees with the set dD up to lower-dimensional pieces. 

5.3. Orientation of boundary. — We now introduce an outer orientation on to 
achieve the stronger property that 0(f?) agrees with the boundary dD as a chain (or 
distribution). For that, fix any x G Gri i (RP+^) and consider some (r, t) E x DiV^). 
Arrange the single eigenvalue A (r) of r and the eigenvalues Ai , A2, . . . , of f as a 

decreasing sequence. (This ordering is uniquely defined except for the set of measure 
zero where some eigenvalues are degenerate.) Following the prescription of Sect. 15.11 
draw the space-like eigenvalues of t as full dots and the time-like eigenvalues as empty 
dots. In view of Fig.[2l if T] (r) > then represent A (r) by •o, else by o». This defines a 
motif C7(r, t) for almost every pair (r, t) . For p = q = 2for example, if Ai > A2 > A (r), 
Ai is space-like, A2 is time-like, and 7](r) > 0, then o{rj) = •o»o. 

The point ^(x;rj) is a point in the boundary of C dD. According to the 

alternating sign of the sum dD = J^sgn(a)(9Dfj the outer orientation of dD at ^{x;rj) 
is directed from the inside to the outside if sgn(a(r,?)) = +1 and from the outside to 
the inside if sgn(a(r,?)) = — 1. Consulting again Fig. [2] we assign to the boundary 
point ^ (x ; r, t) the transverse vector 



where {^*)x-r,t nieans the differential = of the map at the point {x;rj). The 
motif o{rj) switches sign as r}{r) passes through zero; therefore our assignment is 
smooth along the 77 -axis of E^, although discontinuities may occur when the eigen- 
value A(r) hits eigenvalues of t or when eigenvalues of t of opposite type collide. 

In the case of the other component. Fx x D{Vj^), we proceed in the same way, albeit 
with the two boundary coordinate functions ^ and 77 exchanging roles. Thus we assign 
to a point ^ (x ; r, t) for (r, t) e F^x D{Vj-) the transverse vector 



which is again smooth as a function of ^ . 

Now while the distinction between ^ and 77 (or E^ and Fx) is only defined locally in 
X, the outcome of our discussion is free of the ambiguity of exchanging ^ ^ rj. Indeed, 
a glance at Fig. |2] shows that the picture is invariant under the reflection ^ rj) at the 
axis dividing the first quadrant. Hence our procedure gives a globally defined piecewise 
smooth transverse vector field and hence a piecewise smooth outer orientation of the 
boundary ^{S'). Adopting this outer orientation we have 



sign(r](r))sgn(a(r,?))((^,)(^.,^,)(a^) , 



sign(^(r))sgn((j(r,0)(<|)*)(x;r,o(^i7) ^ 



H^)-Y,^^gri{a)dD„ 



as an equality between chains, i.e., in the weak sense: 




(5.3) 
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for any bounded twisted differential form d(0 of top degree. This result holds true, in 
particular, for the integrand CO = iGA,e \dR\ of equation (|3.4I) with the Gaussian 

GA,£(i?)=e-TK^+'^)';f,(i? + iA) 
regularized by the cutoff function Xei^) = e^J^'^^-^^^^")^. 

5.4. Pulling back to S'. — In the next step the left-hand side of the integral (15.31) for 
CO = iGA,e It:(^a) I^^I puUed back to S' C . This results in an iterated integral where 
the outer integral is over the compact manifold Gri i (R^^+^) 3 x and the inner integral 
is over {E^ U F^) x D{V^). It will now be shown that the inner integral over U is of 
exactly the same form as in the case of p = ^ = 1 and hence, for the reasons explained 
in Sect. 14.31 vanishes in the limit e ^ 0. We will do this by inspection of ^*co. 

First we compute the pull back \dR\). Denoting by (^^^)* = (0*)^^ = 

the differential of the inverse map we have 

r{l,(A)\dR\) = l^-\,i^A))^*\dR\. 

With this relation in mind we now compute (^*\dR\. For this we notice that the map 
^ : ^ — > Symg(M'^+^), {x\r,t) ^ rllTc + tild — Ylx) schematically has the variation 

5<l) = {5r)U^ + (50 (Id - n,) + (r - t)5U^ . 

The first two terms contribute factors of unity to the Jacobian when one pulls back the 
translation-invariant measure \dR\ of Symg(]RP+^). 

Now the third term is purely off-diagonal w.r.t. the decomposition W^'^ = K ® ^x' '■ 

5n, e (Hom(y,,y/) ©Hom(y/,yj) nSymg(R^^+0 , 

and may be expressed as a commutator dUx— [oc, 11;^] with a generator a E Lie(Op,q) . 
The two components of 511;,; in Hom(y^,y^^) and Hom(y^^,y^) are related to each 
other by the condition of 5- symmetry. Identifying Vx with its dual V* by means of 
the non-degenerate form 5|\/^x\/i we have an isomorphism Hom(y^,y,^) ^ 
By this isomorphism we may view the multiplication operator dU^ ^ {^ — t) 5Tlx as a 
linear transformation r®lA — lA®t onV^® V^. We thus arrive at the formula 

^*.,.^ jJdR\=Jxirx, tx) dx dr^ dt^ , 

with the Jacobian of the transformation being 

7x(«, b) = |Det(a ldy± - Idy^ 0b) \ . 

Here dx denotes a suitably normalized SOp.g-invariant measure for the compact man- 
ifold Gri i(M^+'^), while drx and dt^ are translation-invariant measures for Symg(yx) 
and D{Vj^). As an observation to be used presently, we note that the Jacobian 

Ua,b) = 'fl I (A (a) -Xk{b)f-^ (a) 77 (a) | 

k=l 

becomes independent of ^ and rj for a E E^U Fx. 
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We now have to contract the density ^*\dR\ with the vector field 0^^(t(A)) and 
make the restriction to S'. Since =0 on and dr] = on , the contraction after 
restriction depends only on the t,- and T]-components of ^^^^(t(A)): 

While these components vary with x , 

A^{x)=B{e-Ae-) , A'^W =5(e+ , 

they are constant along E^^ and F^. Hence, recalling the expression dr^ = dXdt,dr\ 
and assembling factors, we see that {i^f^j^-^\dR\) for fixed x is proportional to the 
form rkdXdX] (with constant coefficient) along E^ and ±LdXdt, (still with constant 
coefficient) along . In other words, the inner dependence (along E^ and F^) of the 
integration form (^*{i^^j^-^\dR\) is the same as in Sect. 14.31 

5.5. Gaussian integrand on boundary. — Our final step is to compute the restric- 
tion to (S' of the pull back (^*Ga,£ of the Gaussian Ga.e ■ For this we decompose 

Tr{R + \Af = Tri?2 + 2i TrAR - TrA^ =: /2 {R) + fi (R) + fo . 

Fixing any x E Gri^i(M^+^) and inserting R = ^{x;r,t) = rllx + t{ld — llx) we obtain 

{rf2){x;r,t) = Tryy)+Try^^{t^) , 

(07i)(x;r,O = 2iTr^,^(An,r) + 2iTr^x(A(Id-n,)0 . 

Next we fix ? G D{V^) and restrict the range of the variable r E Symg{Vx) to ExUF^. 
With this restriction we get 

( * /2) (-^ ; ■ , 1 £^UF, = 2^ ^ + const , 

{(l)*fi){x;;t)\E^ = 2iXTTvM)+2'^riB{e+,Ae+) +const , 

i(l)*fi)ix;-j)\^^ = 2iATrv^^(A) +2i^5(e;,Ae;) +const . 

The coefficients B{ef,Aef) never vanish. Indeed, from A = Y.a=i , ■) we have 

Ae±) = £5((p„,e±)2>0, 

a=l 

by the positivity assumption As > and the non-degeneracy of B. As a consequence 
of the fact that B{ef,Aef) ^ 0, it will suffice for our purpose of taking the limit e ^ 
to examine the cutoff function R t-^ XeiR + '^) at A = . 

By the relations se^ = and se^ = the operator s commutes with the projector 
Ilx = e^5(e^,-) + e^5(e+,-). Therefore the pull back ^*X£ of the cutoff function 
R Q-jTiisR-Rs) separates into two factors, one for and each. We again keep 
t E D{Vj^) fixed and investigate the dependence on r G Symg(y,^). Using the easily 
verified relation H^isr — rs)Ilx = (r) — rj (r))n^ we obtain the expression 
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We finally conclude that 0*(U vanishes for CO = iGA,£ l^^^-j in the limit e — > 0. 
Indeed, we now see from the above expressions that the inner integral over Ej U is 

\Jm. Jr 
and this always goes to zero when the regularization e is sent to zero. 
This shows that ^I{A + tA) \t^Q = and completes the proof of Thm.fT] 
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